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MATHEMATICAL  ANALYSIS  OF  CIRCULAR  CORROSION  CELLS 
HAVING  UNEQUAL  POLARIZATION  PARAMETERS 


INTRODUCnoN 

In  many  corrosion  reactions  the  anode  and  cathode  are  spatially  localized.  This  may 
occur  on  different  surfaces,  as  in  the  galvanic  corrosion  of  dissimilar  metals,  or  on  dif- 
ferent parts  of  the  same  surface,  as  with  localized  geometries  such  as  crevices.  Moreover 
the  corroding  system  often  has  a coplanar  concentric  circular  geometry.  As  pointed  out 
earlier  [1,  2],  examples  include  some  instances  of  pitting  [3J,  crevice  corrosion  under  O- 
rings  or  washers  [4],  and  corrosion  under  barnacles  [5],  under  tubercules  of  corrosion 
products  [6],  or  under  du.st  particles  in  condensed  moisture  films  [7]. 

In  all  such  cases  there  is  a potential  difference  between  the  central  anode  and  the 
disk-shaped  cathode  surrounding  it.  This  potential  difference  may  arise  from  heter- 
ogeneities in  the  solid  phase  (such  as  dissimihu-  metals,  inclusions  in  a base  metal,  or 
discontinuities  in  protective  films)  or  from  heterogeneities  in  the  liquid  phase  (such  as 
differential  aeration  as  in  crevices).  TTius  there  is  a distribution  of  both  electrode  poten-  i 

tial  and  local  current  density  as  one  moves  radially  from  the  center  of  the  anode  out 
toward  the  far  edge  of  the  cathode. 

Gal-Or,  Raz,  and  Yahalom  [8]  have  mathematically  treated  systems  of  coplanar 
concentric  circular  corrosion  cells.  These  authors  analyzed  the  effect  of  various  system 
parameters  on  the  total  current,  and  more  recently  McCafferty  (1,  2]  has  evaluated 
the  distribution  of  potential  and  current  across  such  cells.  'Ihese  treatments  essentially 
extended  to  cylindrical  geometries  the  model  developed  by  Waber  and  coworkers  (9-12) 
in  a series  of  publications  treating  semi-infinite  parallel  electrodes. 

Two  central  features  in  the  Waber  model  are  that  the  ancxle  and  cathode  obey  linear 
polarization  kinetics  over  an  extended  potential  range  and  that  the  anodic  and  cathodic 
slopes  are  equal.  Whereas  the  first  assumption  often  holds  in  expieriments,  the  second 
assumption  rarely  holds,  because  the  anode  is  generally  far  less  polarizable  than  the 
cathode. 

Tbe  case  of  unequal  anodic  and  cathodic  linea  ition  has  been  solved  recently 

by  Kennard  and  Waber  [13]  for  semi-infinite  strips  el  electrodes  under  bulk 

electrolyte. 

This  report  extends  the  Waber  model  of  linear  corrosion  kinetics  to  circular  systems  i 

with  unequal  polarization  parameters.  Equations  are  derived  for  potential  and  current  i 

distributions  and  for  the  total  anodic  current,  and  „jneralized  calculations  are  made.  ■ 

Comparisons  with  experimental  results  will  lie  made  elsewhere.  i 


Manuscript  submitted  February  7,  1977. 
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DESCRIPTION  OF  THE  MODEL 
Hie  Corrosion  Cell 

Hie  cell  geometry  is  shown  in  Fig.  1.  The  anode  and  cathode  outer  edges  are 
coplanar  concentric  circles  of  radii  a and  c respectively.  Hie  electrolyte  thickneu  b is 
aUowed  to  approach  infinity  for  bulk  electrolyte. 


Linear  Polarization 

Hie  cell  potentials  are  shown  in  Fig.  2a,  and  stylized  polarization  curves  are  shown 
in  Figs.  2b  and  2c.  Following  Wagner  [14]  and  Waber  [9-13],  an  important  feature  of 
the  model  is  that  the  polarization  curves  are  linear  in  £ vs  i over  an  extended  range.  As 
pointed  out  by  Kennard  and  Waber  [13],  if  the  plots  are  linear  over  only  a portion  of 
the  curve,  tangent  approximations  can  be  drawn.  Hiiu  the  open-circuit  potentials  and 
are  replaced  by  the  intersections  E^'  and  E^'  respectively  of  the  tan^t  lines  with  the 
potential  axis,  as  shown  in  Fig.  2c. 
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Fig.  2a— Electrode  potentials  across  the  cell. 
and  refer  to  the  open-circuit  potentials  of 
the  anode  and  cathode  respectively. 
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Fig.  2b— Ideal  linear  polarization  curves  over  an 
extended  potential  range 
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Fig.  2c — Linear  approximations  to  the  polariza- 
tion curves.  The  extrapolated  values  and 
replace  and  Ef  respectively. 
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The  linearized  polarization  curves  are  characterized  by  the  Wagner  polarization 
parameters 

^0  = 

for  the  anode  and 

= 

for  the  cathode,  where  a is  the  electrolyte  ^ 
the  dimensions  of  length  (cm),  and  ^ 

The  assumption  of  linear  polarization  oyer  an  extended  potential  range  has  been 
observed  to  be  a reasonable  approximation  in  a number  of  instances.  For  example,  steels 
in  aerated  neutral  to  basic  solutions,  with  or  without  chloride,  displayed  both  anodic  and 
cathodic  plots  which  were  approximately  linear  over  an  extended  range  [15].  Additional 
examples  include  the  behavior  of  copper/steel  couples  in  distilled  water  [16],  the  cor- 
rosion of  tin  in  citrate  solutions  [17],  and  the  corrosion  of  bare  and  coated  aluminum  in 
chloride  solutions  [18].  Other  examples  involve  specialized  geometries,  such  as  the 
pitting  of  aluminum  [19],  or  specialized  conditions,  such  as  the  dissolution  of  mild  steel 
at  high  anodic  overpotentials  in  concentrated  electrolytes  [20] . On  the  cathodic  side  the 
reduction  of  oxygen  on  nickel  in  dilute  H2SO4  [21]  and  of  silver  in  KOH  [22]  display 
linear  regions.  Polarization  curves  for  a variety  of  metals  in  thin-layer  electrolytes  [23] 
display  linear  regions  over  at  least  part  of  the  potential  ranges  for  both  anodic  and 
cathodic  processes. 

In  some  cases  the  linearity  may  be  attributed  predominantly  to  resistance  polariza- 
tion, caused  either  by  iR  drops  through  the  solution  or  by  ohmic  films  on  the  electrode 
surface.  As  pointed  out  by  Stem  and  Geary  [24],  however,  sometimes  the  combined 
effects  of  concentration  polarization  plus  ohmic  drops  interfere  with  activation  polariza- 
tion processes  so  that  a very  short  Tafel  region  is  observed.  Such  cases  often  give 
straight-line  segments  in  E vs  i. 

At  this  point  it  should  be  clear  that  the  model  invokes  linearity  over  an  extended 
potential  range  and  not  merely  in  the  pre-Tafel  region  near  the  corrosion  potential,  where 
the  usual  Stem  and  Geary  [24]  linear  relation  is  valid. 

MATHEMATICAL  ANALYSIS  FOR  BULK  ELECTROLYTE 

The  electrostatic  potential  P(x,  y,  z)  is  given  by  Laplace’s  equation 

V^P(x,  y,  z)  = 0,  (3) 

provided  that  there  are  no  concentration  gradients  in  the  solution,  the  solution  is  electro- 
neutral, and  there  are  no  sources  or  sinks  of  ions  in  the  electrolyte  [25] . 

With  the  circular  geometry  it  is  convenient  to  rewrite  Eq.  (3)  in  cyclindrical  coordi- 
nates using  the  usual  transformations  x = r cos  0 and  y = r sin  0 . The  result  is 


\dE 

I di 


(1) 


dc 


(2) 


conductivity.  The  parameters  and  Eg  have 
, in  the  present  treatment. 
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d^P(r,z)  1 dP(r,  z)  d^P(r,  z)  ^ 

+ = 0,  (4) 

3r2  '•  a--  ^^2 

where  the  potential  P is  independent  of  the  angle  0.  The  general  approach  is  to  solve  for 
P(r,  z)  in  Eq.  (4)  subject  to  appropriate  boundary  conditions  and  then  to  evaluate  the 
local  current  density  i(r,  0)  from  Ohm’s  law  for  electrolytes; 

, (5) 

j = 0 

where  a is  the  electrolyte  conductivity. 


Boundary  Conditions 

The  boundary  conditions  have  been  discussed  in  some  detail  in  a previous  report 
[1] . In  brief  there  is  no  current  flow  across  the  symmetry  line  r = 0,  nor  across  the  cathode 
outer  boundary  r = c.  Thus 

9P(r,  z) 

_ 3r  _ 

and 

^P{r,  z) 
dr 

Also,  the  potential  must  be  bounded  at  the  upper  physical  boundar>'  of  the  electrolyte, 
so  that 

lim  P(r,  z)  < M,  (8) 

2-^00 

where  M is  some  finite  number. 

TTie  general  solution  to  Eq.  (4)  subject  to  the  boundary  conditions  of  Eqs.  (6) 
throu^  (8)  is  [1,  2,  8] 

oo 

P(r,2)  = Co  + ^ (9^ 

where  Cq  and  C„  are  coefficients  to  be  evaluated  later,  Jq  is  the  Bessel  function  of  order 
zero,  and  X„  = x„!c,  in  which  the  are  the  zeros  of  J^(x)  = 0. 


= 0. 


(71 


i(r,  0)  = - a 


dP(r,  z) 


dz 
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Linear  Polarization 


The  remaining  boundary  condition  relates  the  electrode  potential  E{r,  0)  along  the 
metal  surface  vs  some  standard  reference  electrode  to  the  electrostatic  potential  P(r,  0) 
within  the  electrolyte  but  “just  outside”  [26]  the  electrode  surface.  If  and  E^  are 
the  potentials  of  the  polarized  anode  and  cathode  respectively,  at  any  current  density, 
then 


V - P(r,0)  = E{r,0),  (10) 

where  V'  is  a constant  which  includes  the  various  differences  in  electrostatic  potential 
across  the  extra  interfaces  introduced  in  the  measurement  of  a potential  difference  across 
the  metal /solution  interface  of  interest  [27].  Equation  (10)  is  developed  in  Appendix 
A. 


For  the  anodic  branch  in  Fig.  2b, 


slope 


i(r,  0)  - 0 _ 


£(r,  0)  - 

which  after  rearranging  becomes 

E{r,  0)  = + i(r,  0) 

Substitution  of  Eq.  (5)  in  Eq.  (12)  gives 

E(r,  0)  = E°  - a 


di 


dE 


dE 


di 


dE 

3P(r,  z) 

di 

a 

dz 

z = 0 


Use  of  the  Wagner  polarization  parameter  as  defined  in  Eq.  (1)  gives: 

\bP(r,  z)] 


E(r,  0)  = - JC, 


-'2  = 0 


Substitution  of  Eq.  (14)  in  Eq.  (10)  gives 

rdP(r,  2)1 


P(r,  0)  - 


dz 


= V'  - 0 < r < a. 


-J2  = 0 


A similar  expression  holds  for  the  cathode: 

r9P(r,2)1 


P(r,  0)  - £, 


dz 


= V'  - E°  a < r < c. 


2 = 0 


(11) 


(12) 


(13) 


(14) 


(15a) 


(15b) 
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Equations  (15a)  and  (15b)  are  the  final  boundary  conditions  required.  The  indeterminate 
constant  V'  will  vanish  in  the  final  forms  of  the  expressions  to  be  derived  but  will  be 
carried  along  for  mathematical  completeness. 


Evaluation  of  the  Coefficients  C„ 

The  boundary  conditions  in  Eqs.  (15a)  and  (15b)  are  used  to  determine  the  coeffi- 
cients C„  appearing  in  Eq.  (9).  Tlie  approach  is  to  evaluate  Eqs.  (15a)  and  (15b)  using 
the  general  expression  for  P(r,  z)  and  then  to  solve  the  two  simultaneous  equations.  The 
reader  wishing  to  avoid  the  mathematical  details  can  skip  to  Eq.  (35). 

Tlie  general  expression  for  P(r,  z)  was  given  earlier  by  Eq.  (9).  Use  of  Eq.  (9)  in 
(15a)  gives 

OO 

Co  ♦ 2]  * ^aK)Jo(Kr)  = V - E°.  0 < r < a.  (16) 

fl»  1 

If  this  equation  is  multiplied  through  by  rJQ(\^r)  and  integrated  over  the  domain  of 
applicability  (from  r = 0tor  = a),  then 

«'r»0  *'>■•0  n«l 

= (V'-E^)  f rJQ{X„r)dr.  (17) 

•'r-0 

Tlie  first  and  third  integrals  can  be  evaluated  from  a standard  recursion  formula  for 
Bessel  functions  [28] ; that  is 

^jxJl(x)j  = xJq(x),  (18) 

which,  upon  appropriate  variable  change  and  integration,  gives 

JrjQ(Xr)dr  = -^rJi(Xr).  (19) 

TTie  second  integral  to  be  called  /j,  is 

^2  * I Z C„(l*£,X„yJo(\.r)>/o(Kr)dr  (20) 

”^'■•0  n»l 


or 
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^2  ~ 


^ C„(l  + rJQ{\„r)Jfj[X,„r)dr. 

n = I '^r=0 


(21) 


The  summation  can  b<!  split  into  two  cases:  n - m,  and  n =#  m.  Thus 


^2  = f 

'^r=0 

+ I r^(X„r)7o(X„r)dr. 

«=1  •'r=0 


(22) 


Substitution  of  Eqs.  (19)  and  (22)  back  into  Eq.  (17)  thus  yields 
1 

Co  T-  0'^l(>^m0)  + +,C„X„, ) rj|(X^r)dr 

m -/r=0 


X!  C„(l+XuX„)  rJo(X„r)^(X,„r)c(r 

n=l  -’r=0 


n ^ rn 


= (V'-E")^o.y,(X„,a). 

The  second  intejjral  in  Eq.  (23)  is  a Lommel  integral  |29]  : 

fjc{/3c/„(cvjr)Jn(|j;c)  - aJ,; (o;Ar)J„(i3;c)} 


(23) 


(0!^-  /j2)  f xJ„(oac)J„{iix)dx  = 
1 


Xo 


(24) 


where  the  primes  denote  differentiation  with  respect  to  the  whole  aqjument  and  not  just 
X.  Thus,  when  n m, 

I rXo(^„d</o(X^r)dr  = |X„J,  (\.«l-^o(^m“»  " (^n,“)l  ■ 

•^n  X“  - X^ 


(25) 


When  n = m,  the  mU-Kral  becomes  frJ^(X„,r)  dr,  which  is  the  remaining  inU'gral  to  U’ 
evaluated  in  Eq.  (23).  However,  when  n = m,  the  right-hand  side  of  Erj.  (25)  gives  0/0, 
so  that  rHospital’s  rule  must  be  used.  In  this  ca.se  the  ruimerator  luid  denominator  are 
differentiated  with  ivsfK'ct  to  X„,  and  then  X„  is  allowed  to  approach  X^ . The  n>suU 
(omitting  several  steps)  is 
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f <^4 

*'r=0 

Use  of  Eqs.  (25)  and  (26)  in  Eq.  (23)  gives 


(26) 


1 

CqT  aJj(X^a)  + + .CgX^)  (</Q(X^a)  + Jj(X^a) 


+ ^ C„(l+£aX„)__ [X„Ji(X„a)Jo(X^a)  - X^Jo(X„a)J](X„a)] 

n=l 

n^m 


= (U'-fO)  Ajj(X„a). 


(27) 


Equation  (27)  is  one  of  the  two  simultaneous  equations  to  be  solved  for  the  set  C„. 
The  second  equation  follows  from  the  boundary  condition  on  the  cathode  given  in  Eq. 
(15b).  The  approach  is  the  same  as  has  just  been  completed.  P(r,  0)  and  dP{r,  z)ldz  at 
2 = 0 are  evaluated  from  Eq.  (9),  so  that  Eq.  (15b)  becomes 


Co  + XI  C„(l  + £,X„)Jo(^n'-)  = V - E^,  a < r < c. 


(28) 


1 


Again  the  equation  is  multiplied  by  rJoi^m'')  dr  and  integrated  over  the  appropriate 
limits,  which  in  this  case  are  from  r = a to  r = c: 

C C 

Co  r rJo(-K^r)dr  + f X] 

*^r=a  '^r=a  n^l 

= (V'-E°)  r rJo(Xrr,r)dr. 

^1"  3:  /» 


(29) 


The  first  and  third  integrals  can  be  evaluated  using  Eq.  (19): 

C °°  ~c 

— ^ [cJi(X^c)  - a Ji(X,„a)]  + Y[  C„(l+£cX„)  | rJQ(XnT)jQ{\„ 

„=1  -'r-Q 


r)  dr 


V - E" 


[cJj(X^c)  -aJi(X,„a)]. 


(30) 
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By  definition  X^c  •;=  and  = 0.  Again  the  summation  can  be  split  into  two 

cases,  so  that  Eq.  (30)  becomes 

c 

- + JC,X^)  rj|(X„r)dr 

•'r-a 

« -c 

n-l  *^'•■0 


(V'-E°) 


aJi(Xma). 


The  integrals  in  Eq.  (31)  can  be  evaluated  as  before  from  Eq.  (24).  A simpler  approach 
is  to  add  and  subtract  /“  rJ^(X„r)  dr  to  the  second  term  and  to  add  and  subtract 
/“  rJo(X„r)Jo(X^r)  dr  within  the  summation  sign;  Then  Eq.  (31)  becomes 

--^aJl(X^a)  + C^(l  + £c>^;„)  f f rJ^(\„r)dr 

|_-'r-0  *'r«0 

r r‘  r“ 

+ y C„(l+XcX„)  rJQ(\^r)JQ{\^r)dr-  rJQ(\r)jQ(\^r)  dr 

;|ri  L-^r-o  Vo  j 


The  integrals  involving  the  entire  interval  from  r = 0 to  r = c are  the  usual  orthogonality 
relations  [30], 


p rdo(X„r)do(X„r)dr  = 


0 ,n¥=m. 


M (33) 

[«^o(Vc)l  ,n  = m. 


and  the  integrals  from  r = 0 to  r = a have  already  been  evaluated  per  Eqs.  (25)  and  (26), 
so  that  Eq.  (32)  reduces  to 
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c r 2 2 1 

- a</](X„,a)  + C,„ ( 1 + C(.X„, ) |— </q  (X,„c)  ^ [</q (X^a)  + Jj  (X^a)]| 


+ y C„(l  + .CpX^  ) J — 


\1  - 


r aJi(X,„a). 


Equations  (34)  and  (27)  are  thus  two  simultaneous  equations  in  Cq  and  C„. 
Addition  of  the  two  equations  eliminates  Cg-  The  result,  after  considerable  algebra,  is 


1 + 


' ' ' _J  c . n = l 


fa 


7)’ 


where  and  are  defined  by 


_ - 


Vnm  = ^ ^ '^1  p cfO  cj  ■ X„  cj  1 V""  C, 

X n ~ ^ m 


Equation  (35)  thus  generates  a series  of  equations,  say  k of  them,  where  m is  fixed  in  turn 
from  1 through  k.  These  k equations  are  solved  simultaneously  to  give  the  coefficients  C„ 
from  n = 1 to  k. 

The  indeterminate  constant  V cancels  out  in  the  generation  of  the  set  of  C„. 
^\lso,  when  = £,  Eqs.  (35a)  and  (35b)  reduce  to  the  previous  case  [1,  2] ; 


l°\t  ! “ ''i 


c„  


where  the  dummy  variable  m has  been  replaced  by  the  more  general  n. 
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Evdlualion  of  the  Coefficient  (^o 

'I’he  remaining  coefficient  Cq  can  be  most  conveniently  evaluated  by  going  back  to 
the  original  pair  of  equations;  Eqs.  (16)  and  (28).  The  approach  is  to  multiply  through 
in  both  equations  by  rdr  and  then  to  integrate  over  the  appropriate  limits.  The  opera- 
tions are  straightforward,  so  that  it  is  not  necessary  to  detail  the  proof  here.  The 
resulting  two  equations  are 


(37) 


n*  1 


and 


c2 


2 0 o®  / 

- 2]  C„(1+£,XJ^  Ji(X„a)  = (V' -£,")( 

^ M=1  ” 

Adding  Eqs.  (37)  and  (38)  and  solving  for  Cq  gives 

Co  = - 2-(£a-£c)t 


(38) 


n = l 


(39) 


Electrode  Potential 

The  electrostatic  potential  P(r,  z)  at  any  point  in  the  electrolyte  is  given  by  Eq.  (9): 


oo 

L/  r\  -Xn^/c 
^n’^oi^n~p 
n=l  ' ' 


(9) 


Use  of  Eq.  (39)  for  Cq  in  the  above  gives 


P(r,z)  = V' 


2a 

2 


n = l 


OO 

“V 

n = 1 


z/c 


(40) 
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The  relationship  between  the  electrostatic  potential  P(r,  0)  in  the  electrolyte  near 
the  metal  surface  and  the  electrode  potential  £(r,  0)  vs  a standard  reference  electrode 
is  given  by  Eq.  (10).  Use  of  Eq.  (40)  with  2 = 0 in  Eq.  (10)  gives 

/ \ 

E{r,  0)  = ^ Vc  ^ Z]  (■''"  c ) ' 


(41) 

where  the  C„  are  determined  from  Eq.  (35).  Again  the  indeterminate  constant  U' 
vanishes  in  the  final  expression. 

Local  Current  Density 

The  local  current  density  i(r,  0)  along  the  metal  surface  is  relab^d  to  the  electrostatic 
potential  P(r,  z)  by  Eq.  (5).  Performing  the  differentiation  on  P(r,  2)  as  given  in  Eq.  (40) 
and  inserting  the  result  in  Eq.  (5)  yields 


i(r,  0)  _ 1 V 


- Zi  ^n^n'^oi^n  )- 


n = l 


(42) 


where  again  the  set  of  C„  is  determined  from  Eq.  (35). 


Total  Anodic  Current 

The  total  anodic  current  is  related  to  the  local  current  density  by  |1,  2,  8] : 

,a  ,-271 

J 

■r=0  *^0  = 0 


-a  ,-277 

Wc  = I I Hr,0)rdrd0 


(43) 


or 


^anodic  = 277  / .(r,0)rdr 

*'r=0 


(44) 


Use  of  Eq.  (42)  in  (44)  gives 


I 


•'r-O  ^ ' 


anodic 


— ) r dr 


(45) 


or 


I 


anodic 


n-1  "7“0 


■ dr. 


(46) 
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The  integral  can  be  evaluated  from  Eq.  (19),  so  that  the  result,  omitting  a few  steps,  is 


where  again  the  set  of  C„  is  determined  from  Eq.  (35). 


MATHEMATICAL  ANALYSIS  FOR  A THIN-LAYER  ELECTROLYTE 

If  the  electrolyte  is  a thin  layer  of  height  b instead  of  bulk  liquid,  the  boundary 
condition  given  by  Eq.  (8)  is  replaced  by  the  requirement  that  there  is  no  current  flow 
across  the  outer  boundarj'  of  the  electrolyte: 


9P(r,  z) 
3z 


The  other  boundary  conditions  are  the  same  as  for  the  bulk  case.  The  general  solution 
of  Laplace’s  equation  subject  to  the  restrictions  of  Eqs.  (6),  (7),  and  (48)  is 

P(r,  z)  = Cq  + ^ C„  cosh  -^(9-2)  (49) 

The  coefficients  Cq  and  C„  are  evaluated  from  Eqs.  (15a)  and  (15b),  as  was  done  for  the 
case  of  bulk  electrolyte.  The  procedure  is  exactly  the  same  as  for  the  bulk  case;  hence 
only  the  results  are  listed  below. 

The  coefficients  C„  are  given  by  the  systems  of  simultaneous  equations 

J -I  r i 1.  / '^o('^m)  ^mfa\^  I . 1 

4 '(•c  p 4anh  — j ^ ~ ^7/  f 


\2  I l'”7/  J b\ 

co.h  cosh  J 


Z 7 Hi  r">  7 ’ 


X 


where  again  and  are  defined  by  Eqs.  (35a)  and  (35b)  respectively  and  m takes 
on  the  values  1 through  k successively. 


The  constant  Cq  is  given  by 
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c-  = V"  - - E, 


< t<;  - C,)  Cn  sinh  i^x„ 


f;(r,  0)  = Ei;  + Ja;'.'  + — (t„  - L,|  ^ C„  smh  {x^  ^ 


- }]  C„  cosh  ix„-]jo 


The  local  current  density  is  given  by 


i(r,  0)  1 


= r sinh  a:„ -bo  •*« 


and  the  total  anodic  current  is 


I ^ 


When  b-*°o,  C„  (thin-layer)  X sinh  {x„blc)  ->■  C„  (hulk),  so  that  the  expressions  for  thin- 
layer  electrolyte  reduce  to  the  corresponding  equations  for  a bulk  electrolyte  for  large 
b. 

The  various  expressions  for  the  thin-layer  and  bulk  cases  are  summarized  in  Table  1. 


PREVIOUS  CASE  OF  EQUAL  POLARIZATION  PARAMETERS 

When  = £,  the  preceding  results  for  bulk  and  thin-layer  electrolytes  reduce 


E(r,0)  = - U"  + 


OO  J 1 ix„  — ^ . . 


E McCAFFERTY 


TabU'  1 Summarv'  of  Relationships  for  Coplanar  Circular  Electrodes  under  a Bulk  and  a 
Thin-Layer  Klec  lrolyU*.  The  hypt^rbolic  functions  m the  upper  part  of  the  brackets  apply  to 
thin  layers,  and  an*  n*placed  by  1 for  a bulk  electrolyte.  The  numbers  in  parentheses  refer 
to  equation  numliers  in  the  U*xt. 

(Wffich*nts  C„ 


Electrode  Poh*ntial 
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and'. 


i(r.  0) 
o 


CO 


where  in  each  of  the  above 


bulk  e'  •'  trolyte, 
thin  layer. 


(oGi 


(57) 


(58) 


These  equations  are  the  same  as  those  previously  reported  except  that  the  signs  of 
electrode  potentials  in  Eqs.  (55)  through  (57)  now  conform  to  the  convention  that  the 
noble  direction  is  the  more  positive. 


NUMERICAL  EVALUATION 
Values  of  and 

In  general  the  anode  and  cathode  have  different  polarizabilities  (the  two  electrode 
potentials  respond  differently  to  the  passage  of  current).  In  many  instances  the  anode  is 
the  less  pol^uizable.  This  is  illustrated  by  many  electrode  kinetic  studies  carried  out 
under  carefully  controlled  conditions.  With  iron,  for  example,  in  a variety  of  deaerated 
electrolytes,  anodic  Tafel  slopes  of  30  to  80  mV/decade  have  been  observed,  while  the 
cathodic  slopes  were  120  mV/decade  [31-34],  Other  metals  in  the  iron  group  (nickel 
and  cobalt)  have  been  observed  to  behave  simil2U'ly  [35]. 

To  cite  two  more  examples,  cadmium  undergoes  self-dissolution  to  Cd^^  in  acids  by 
two  consecutive  single-electron  transfer  reactions,  and  indium  goes  to  In*  '^  through  three 
consecutive  single-electron  transfers.  The  observed  anodic  Tafel  slopes  are  40  to  50mV/ 
decade  (0°C)  for  cadmium  [36]  and  22  mV/decade  for  indium  [37],  in  good  corre- 
spondence with  the  theoretical  values  of  2.303  RT/(3/2)F  and  2.303  RT/(5/2)E 
respectively.  The  hydrogen-evolution  reaction  on  both  surfaces  gave  cathodic  Tafel  slopes 
of  115  mV/decade  and  120  mV/decade  respectively,  indicative  of  a single-electron  trans- 
fer characterized  by  a theoretical  value  of  2.303  RT/(1/2)F. 

In  more  practical  situations  where  conditions  are  not  as  well  defined,  Tafel  behavior 
is  not  always  observed,  but  instead  polarization  curves  sometimes  display  segments  which 
are  approximately  linear  in  current  (rather  than  the  logarithm  of  current)  over  a con- 
siderable potential  range,  as  discussed  earlier.  In  these  instances  the  cathode  often  again 
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is  the  more  polarizable.  Figure  3a  shows  a schematic  Evans  diagram  138]  for  a bimetallic- 
couple  under  cathodic  control.  Other  possibilities,  however,  include  anodic  control 
(Fig.  3b)  and  mixed  control  (Fig.  3c).  This  last  case  would  approximate*  earlier  trc'at- 
ments  (1,  2,  8]  for  = £^. 

Values  of  Wagner  linear  polarization  curves  compiled  [1,  2)  from  the  literature  in- 
dicate that  is  generally  of  the  order  of  1 to  10  cm  while  £^  is  usually  10  to  100  cm, 
fdthough  there  cire  exceptions.  In  data  tabulated  by  Gouda  and  Mourad  [15]  for  steel 
in  a variety  of  neutral  to  basic  solutions  both  with  and  without  added  chloride,  the 
cathodic  slope  IdE/dil  varied  from  1.5  times  to  approximately  20  times  the  ajiodic  sloi>e 
but  with  most  ratios  in  the  range  of  5 to  10.  Specific  conductivities  were  listed  for  only 
three  solutions,  for  which  values  of  £^  axe  calculated  to  be  0.9,  1.9,  and  48  cm,  with 
corresponding  £^  values  of  7.5,  12.2,  and  193  cm  respectively. 


TOTAL 

CURRENT 

‘11 

V.  1 
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\, 
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- 
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Ei 

iO.OOl 

(b) 

■-  1 
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£'c  = 0.1  : 

Er  ' E 
11.00) 

E2 

(0  001 
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II 

E," 

E° 

Fig.  3— Schematic  representations  of  cath- 
odic, anodic,  and  mixed  control  with 
Xq  fixed  at  1 cm  in  each  case 
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Hulk  Kit'd  roly  U* 

For  the  ealfiiUitiuiis  in  this  seelion,  the  anodie  polarization  piU’anii'ter  is  fixed  at 
Cy  -■  1 em.  Fixture  1 shows  tht'  potential  distrihution  A(r,  0)  for  Cy  = 1 cm  and  = 10 
cm  for  fixetl  values  of  /•.y*  = 0.00  V and  = 1.00  V.  The  coefficients  C„  were  calculated 
up  to  f'n)Q  usintt  tht'  system  of  simultiuieous  equations  tteneraled  by  Kq.  (35).  'Fhese 
simultiuieoiis  equations  were  solved  usin^j  the  C'lX'  3300  eompuLer,  and  the  coefficients 
were  then  substituted  in  Ft].  (11)  to  obtain  the  electrode  poU'ntial  distribution.  Con- 
vergt'tue  was  ;Ls.ses,sed  by  numerical  evaluation.  Thi'  computer  program  is  given  in 
Appendix  B. 

F’igure  4 also  shows  potential  distribution  plots  for  ,Cy  = = 1 cm  and  ,C„  = .C^  = 

10  cm,  calculatt'd  from  Eq.  (55).  It  is  evident  that  the  potential  behavior  of  the 
electrodes  with  unequal  luuxlic  ;uid  cathodic  polarization  parameters  cannot  be  deduced 
from  the  two  separate  curves  for  the  equal  polarization  piu'ameU'rs. 

The  com'sponding  current  density  curves  for  the'  three  sysU'ms  are  shown  in  Fig.  5. 
It  is  seen  that  the  values  for  the  case  of  unequal  parameters  are  intermediate  between  the 
two  cases  where  ,Cy  = = 1 cm  and  ,Cy  = .C^  = 10  cm. 
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Fig.  4 • Comparison  eleci.'ojH  p jp-niial  distributions  for 
equal  and  unequal  polari/.ation  parameters  with  bulk 
electrolyte  (anode  radius  a - 0 5 ern,  cathode  radius  r * 1 .0 
cm  - 0.00  V,  and  ■ 1 00  t'l 
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0 8 • 

ANODf  CATHODE 


RADIUS,  r (cm) 


Fig  5— Current  distributions  corresponding  to  the 
electrode  potential  distributions  in  Fig.  -1 


I 

I 


Figure  6 shows  the  potential  distribution  calculated  from  Eq.  (41)  with  n = 100 
for  a fixed  value  of  = 1 cm  but  with  variable  £^.  Corresponding  current  distribution 
curves  calculated  from  Eq.  (42)  an*  shown  in  Fig.  7.  When  = 1 cm  and  X^.  = 0.1  cm, 
the  galvanic  couple  is  under  anodic  control,  as  depicU'd  in  Fig.  3b,  and  the  electrode 
potentials  across  the  metal  surface  of  both  components  are  polarized  up  near  the  potential 
of  the  uncoupled  cathode. 

When  X^  » X^,  such  as  = 1 cm  and  X^  = 100  cm.  the  system  is  under  cathodic 
control,  as  illustrated  schematically  in  Fig.  3a.  For  this  case.  Fig.  3a  predicts  that  the 
electrode  potential  would  approach  the  values  of  the  o(ien-circuit  potential  of  the  anode 
and  the  current  would  be  much  less  than  for  the  case  of  amxiic  control  (for  fixed  X(,). 
These  expected  trends  are  verified  in  the  results  of  the  numerical  analysis  shown  in  Figs. 

6 and  7. 
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0 ' ’ ' 0^5  ’ ' ‘ * 10 


RADIUS,  r Icml 

Fig.  6— Distribution  of  electrode  potential  across  circular 
cells  under  bulk  electrolyte  with  £q  fixed  at  1 cm,  com- 
bined with  various  values  of  (anode  radius  a = 0.5  cm, 
cathode  radius  c ■=  1.0  cm,  = 0.00  V,  and  E°  = 1.00  V’) 


Results  in  Figs.  6 and  7 for  = 1 cm  also  provide  a check  on  the  consistency 

of  the  present  method  with  the  previous  relationships  for  equal  polarization  parameters. 
Both  current  and  potential  distributions  calculated  from  the  set  of  C„  resulting  from  Eq 
(35)  with  Xg  and  X^  both  equal  to  1 cm  agree  with  the  results  obtained  from  Eqs.  (55) 
and  (56). 

One  additional  trend  can  be  seen  in  Figs.  6 and  7.  For  this  system  of  fixed  L^,  the 
more  polarizable  the  cathode  (the  larger  X^.),  the  more  uniform  the  potential  and  current 
distribution. 


The  tot£il  tinodic  current  was  calculated  from  Eq.  (47)  for  the  systems  with  X^  fixed 
at  1 cm  with  variable  X^.  Results  are  listed  in  Table  2. 

The  total  smodic  current  can  also  be  calculated  from  the  schematic  Evans  diagrams 
shown  in  Fig.  3.  For  the  anodic  branch 


u 


E - 

^conr 


(59) 
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Fig.  7— Distribution  of  current  acro.ss  circular  cells  under  bulk 
electrolyte  with  Cg  “I  ' '’t"  combined  with  various 

value.s  of  £<■  (anode  radius  u 0.')  cm,  cathode  radius  c = 1.0  cm, 
/•;"  - 0.00  V,  and  F.';.'  = 1.00  V) 


and  frtr  tho  (•athodir  brani  h 


df; 

d/' 


* con 


(60) 


whfTo  ^ rnrr  oorrosion  rtirrt'nt  (the*  total  antalic  currtMit  rt'fpiTod  to  carlior  as 

I,  ^ ^C  ^ '■(•''V  •''here  and  .-1^.  artt  the  aroa  of  anodo  and  cathode  respectively. 

Us«'  of  K(|s.  ;i)  and  (2)  m Kqs  (59)  and  (69)  ffi-.e.r 


2? 


N'KI,  RKl’OUT  H1U7 


'I'ahlf  2 — C'ompanson  of  Tot;il  Currcnl  t';iJi  ulalfci  from  Kvan>.  l)i;i*'rams  ;ukJ  From  Summa- 
tion of  Curronl  Dislnhution  Curves  for  ('iri'ular  Couples  I'liticr  Bulk  Fleetrolyie  (;uiode  radius 
a = 0.5  cm,  cathode  radius  c = 1.0  cm,  /-.'y  = 0.00  C,  and  = 1.00  V) 


C. 

(cm) 

f-c 

(cm) 

/corr/<C 

CalculaU-d  frcrni  Evans 
Diagrams:  Eep  ( 01  ) 

AmeclK- 

Calculated  from 
E(^^.  ( .15  1 :uid  ( 47  1 

" V'  " 

0.1 

0.700 

0.007 

1 

0.5  89 

0.485 

2 

0.471 

0.401 

10 

0.181 

0.109 

20 

0.103 

0.099 

50 

0.0445 

0.0437 

100 

0.0228 

0.0227 

10 

10 

0.0589 

0.0570 

100 

0.0181 

0.0179 

100 

100 

0.00589 

0.00589 

^corr 

U 


(01  ) 


when*  again  has  the  same  meaning  of  /anodic  ^'0-  ClT). 

Values  of  /anodic/^^  calculated  from  Kep  (47)  lue  also  listed  in  lahk'  2.  I'hese  caku- 
lated  values  agree  with  the  results  from  the  computer  analysis  for  = 1 cm  coupled 
with  cathodic  values  of  = 50  cm  and  = 100  cm,  where  the  current  distribution  is 
uniform,  as  shown  in  Fig.  7.  There  is  disagreement  between  the  n'sults  of  the  Kvans- 
diagram  analysis  and  the  computer  analysis  for  tho.se  systems  where  there  is  a nonuniform 
distribution  of  current,  and  this  divergence  is  greater  the  more  pronounced  the  localized 
attack  at  the  anode/cathode  juncture. 

Results  for  .C^  = C^.  = 10  cm  ;uid  C.^  = 1^.  = 100  cm  an-  also  included  in  Table  2. 
Current  distribution  plots  published  in  ;m  earlier  report  1 1 1 were  nearly  uniform  lor  the 
former  system  and  exactly  so  for  the  latter.  The  current  distribution  for  C,,  = 10  cm  and 
Xj.  = 100  cm  was  calculated  from  Eij.  (12)  and  was  also  observed  to  Ix'  uniform  ( I'ht'  plot 
is  not  shown  here.)  Thus  there  is  good  agreemenl  betweem  Fvans-tiiagram  analyses  and 
computer  calculations  for  the  cases  where  (here  is  a uniform  current  distnbuti<m. 

'rhus  the  classic  Evans  polariz.ation  diagrams  cannot  be  used  to  accurately  predict 
the  value  of  galvanic  currents  unless  the  anode  and  cathode  components  each  bi  have 
uniformly. 
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'ITiin-Layer  Electrolyte 

Figure  8 shows  the  electrode  poU’iitial  A’(r,  0)  for  = 1 cm  and  = 10  cm  for 
different  electrolyte  thicknesses.  The  coefficients  C\  were  computed  uj)  to  from 

Eq.  (50)  and  were  used  in  Eq.  (52)  to  obtain  the  poUmtial  distributions.  The  computer 
program  for  thin  layers  is  given  in  ApjK’ndix  C. 

Figure  8 shows  that  the  potential  distribution  is  almost  uniform  for  bulk  electrolyU-, 
but  with  thin  layers  most  of  the  polarization  takes  place  near  the  anode /cathtxJe  juncture. 
The  anode  center  and  cathode  outer  edge  are  virtually  unaffecU'd  by  the  presence  of  each 
other  for  the  thinnest  electrolyte  of  0.001  cm. 

'Fhe  corresponding  current  distributions  are  shown  in  Fig.  9.  The  local  current 
densities  were  calculated  from  Eqs.  (50)  and  (53)  with  n = 100,  except  near  r = 0.0  and 
r = 0.5,  where  125  terms  were  used.  For  the  thin  layers  there  is  a geometry  effect  in 
which  the  corrosion  attack  is  concentrated  near  the  tinode/cathodc  boundar>\ 
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Fig.  9— Current  distribution  for  = I cm  and 
for  different  electrolyte  thicknesses.  The  cell  parameters 
are  the  same  as  in  Fig.  8.  (Limits  between  which  the  local 
current  densities  oscillate  as  computed  from  Eq.  (53)  are 
indicated  for  the  anodic  points.  Limits  are  not  shown  for 
the  cathode  but  are  approximately  half  the  range  of  tte 
anodic  points) 


Figure  10  shows  the  total  anodic  current  calculated  from  Eq.  (54)  for  two  different 
combinations  of  and  .£^.  In  both  cases  the  total  current  approaches  values  for  the  bulk 
for  electrolyte  thicknesses  of  approximately  0.1  to  0.3  cm. 


Figure  11  compares  the  potential  distribution  for  Xq  = 1 cm  and  = 10  cm  to  the 
cases  of  equal  polarization  parameters:  X^  = X^,  = 1 cm  and  Xq  = X^  ■ 10  cm  for  an 
electrolyte  thickness  of  0.001  cm.  Figure  12  shows  a similar  curve  for  X^  = 10  cm  and 
= 100  cm.  In  each  case  the  potential  distribution  for  the  system  of  unequal  parame- 


ters is  not  related  in  a simple  manner  to  the  individual  distribution  curves  for  each  of  the 
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Fig.  10— Total  anodic  current  computed  as  a function 
of  electrolyte  thickness  for  two  different  combina- 
tions of  and  Lc  (anode  radius  a * 0.5  cm,  cathode 
radius  c = 1.0  cm,  = 0.00  V',  and  = 1.00  V) 


Fig.  1 1 —Comparison  of  electrode  potential  distribution 
for  equal  and  unequal  polarization  parameters  for  a thin- 
layer  electrolyte  of  thickness  b * 0.00 1 cm  (anode  radius 
a = 0.5  cm,  cathode  radius  c * 1.0  cm,  = 0.00  V,  and 
Ey  * 1.00  V') 
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Fig.  12— A second  comparison  of  electrode  potential  dis- 
tribution for  equal  and  unequal  polarization  parameters 
for  a thin-layer  electrolyte  of  thickness  b = 0.001  cm 
(anode  radius  a = O..!)  cm,  cathode  radius  c = 1.0  cm,  E”  = 
0.00  V.  and  £,”  = 1.00  V) 


SUMMARY 

A mathematical  model  has  been  developed  to  describe  the  distribution  of  potential 
and  current  across  circular  corrosion  cells  having  unequal  anodic  and  cathodic  lineai' 
polarization  parameters.  This  analysis  is  applicable  to  systems  of  bimetaJlic  galvanic 
couples  or  to  systems  with  a localized  geometry  effect,  as  in  pitting  corrosion. 

ITie  potential  distribution  in  a system  having  unequal  anodic  and  cathodic  polariza- 
tion parameters  is  not  related  in  a simple  manner  to  the  separate  distribution  curves  for 
the  two  cases  where  the  polarization  parameters  are  t>qual. 

For  bulk  electrolyte  the  value  of  the  electrcxle  potentials  depends  on  whether  the 
system  is  under  anodic,  cathodic,  or  mixed  control.  Current  distribution  is  more  unifomi 
for  the  more  polarizable  combinations  of  electrodes.  Thus  the  total  corrosion  < urrent 
calculated  from  the  Evems  diagram  is  in  error  if  the  individual  current  distributions  are 
not  uniform. 

In  thin-layer  electrolytes  there  is  a geometry  effect  in  which  the  electrode  polanza- 
tion  and  current  flow  is  concentrated  near  the  anode/cathode  juncture.  In  a typical  sys- 
tem the  tendency  toward  bulk  behavior  occurs  at  about  0.1  to  0.3  cm  (1000  to  3000 
pm). 
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Appendix  A 


RELATIONSHIP  BETWEEN  THE  ELECTRODE  POTENTIAL  AND  THE 
ELECTROSTATIC  POTENTIAL  AT  THE  METAL/SOLUTION  INTERFACE 


■Vs  puinU'd  out  by  Bockris  and  Reddy  |A1),  it  is  impossible  to  measure  the  electrode 
potential  of  u metal/solution  interface  without  introducing  additional  extraneous  inter- 
faces during  the  measurement  process.  This  is  illustrated  in  Fig.  Al,  where  the  electrode 
(Kitentiais  of  the  coplanar  anode  8ind  cathode  are  to  be  measurt'd  versus  the  reference 
electrode. 


Fitf.  A 1 — Method  of  measurinj»  the  el»*ctrode  potentials  of  a 
coplanar  anode  (A)  and  cathode  (C) 


In  circuit  1 the  measumd  electrode  potential  of  the  anode  Vo  the  reference 
electrode  is  related  to  the  potential  differences  across  the  various  interfaces  by 

0)1  + [F(r,'o)-0„^f|  + (C-RK -O.w;  » = <A1) 

where  0,^  is  the  electrostatic  potential  “just  inside”  the  metal  A and  P(r,  0)  is  the  electro- 
static potential  in  the  solution  “just  outside”  tlu-  metal  |A2).  Similarly  c!>Rpf  and 
refer  respictively  to  the  electrostatic  potraiti il  just  inside  the  solid-phase'  reference  elec- 
trode and  just  insicle  the  connecting  win'. 
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According  to  Bockris  and  Reddy  (Al|  the  potential  difference  across  a nonpolari- 
zable  interface  such  as  /solution  is  a constant,  so  that 

<!>A  ^('■.0)  *<^Ref  ^ 

where  V'  is  a constant.  Use  of  this  definition  of  V'  in  Eq.  (Al)  gives 

V'  - P{r,0)  = 0 < r < a.  {A 

Similarly  measurement  of  the  electrode  potential  of  the  cathode  in  circuit  2 gives 
(0^-</>c)  + I0c 0)1  + [P(r,O)-0^f]  + (0ref-0.wi>  = f'c- 

The  0c  terms  in  Eq.  (A4)  cancel,  so  that 

V'  - P{r,  0)  = E , a < r c.  (A 


V - P{r,  0) 


£„,  0 < r < a 


E , a < r < c 


Or  denoting  the  electrode  potential  along  the  metal  surface  as  E(r,  0)  gives 


V'  - P(r,  0)  = E(r,  0), 


which  is  Eq.  (10)  in  the  main  text 
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Appendix  B 

COMPUTER  PROGRAM  FOR  COPLANAR  CONCENTRIC  CIRCULAR 
ELECTRODES  WITH  UNEQUAL  POLARIZATION  PARAMETERS 
UNDER  BULK  ELECTROLYTE 


yNtGBULK 

C 

y TMS  Pwf-ciWAM  LUKFyTfcS 

C (1)  CLFBEMT  tIS^NIaL^IE^ 

C (2)  FEIEMIAL  UbT«lfcl,l  ItN 

C U)  TSlAy  AKIiEIC  CuBFfcNT 

C FFF  Ce^lCfc^THlC  CIRCULAR'  fclfcCTCErES  CevERt"  RY  3^1*^  FlEC^RFLYTC 

c fPR  THt  c*bE  'where  anELIC  AKC  tAlRtnc  -aGnEH  Rtiy  a r I 2 a T 1 oN 

c EARARtrfHS  ARE  AtT  ECEAy, 

C 

t AiRAtlub  ai'  aNUCE 

t LaRALIuS  aE  catrehe 

y RsLiaTAAyt  alQiyG  RAiiiLb 

c RRrB/r 

y LAAAAF.Jly  wAGNER  R£laR  iZA't  IFN  FaRaREIEr 

y yCACATHEyly  .agree  peuARI^aiifn  farafetEh 

y *(F)JMTH  2E«a  af  RkSSEL  FLNCTIEA  EE  £HDEH  JIARP 

y F (R,c )a1nTEHFACIAL  R£TEaTIAL  ALEKG  TRE  metal  SJRFACE 

c E (R  , 0 )sFd  TENELsELECTErE  PfIENTJAL  ALEMG  TE^E  METAy  SURFaCE 

y E < R , 0 ) :C0*4bTANT-E  (R  , C ) I rITE  TRE  rtNETAM  CANCELylNr.  AyT  IN  1 

y final  EAEiRESblaN,  sa  CENblANT  in  EEEeCT  ran  BF  set  EaJAL  T£  2FRR. 

y 

y THE  y'-'EF  E lylESTE  CO  ANL  OEaRN  are  LEEInEL  IN  The  EollTIrInG  FLUAYIRn 

c 

C K 

y F(R,C)  S y . bLR  (C  . (x(N)«R/:)) 

y 3 N=i  N c 

y 

C I’ETALsTgrAy  A.,'-UC  CLRFEET  CIVILEL  dY  TrF  CSNOuc  T [ v I Tr 

y ILECALstaGAL  CURRENT  TENSITY  rlvlTEC  HY  THE  CdNDuCT I v I TV 

c 

real  la lyC • 1 tpT AL , Il£C Ay 

LImENSIEN  A(lD0),Y(10C.10C).b(lCC.lC0) 

FlAj,lAl&U26b36 
R = 100 

1 feat  IC.A.C.yA.LC.EA.FL 

13  E PRR  AT  (6E  1,' , 0 ) 

PRINT  11,A,C,LA.lC,EA,EC 

11  ECRhaT  (1H1,SX,2RA1  FIO.S.SX.'RC*  E 1C,S,5»,3"LAS  E yy . S , Sx  , JhlC=  Ei 
10 .b .i* . JmE A=  E 1C .S .bx , 3REC A FiO.t.////) 

PRINT  12 

12  EERmaT  (1X..GEHERAT1EE  EE  Tee  System  EF  S I muL  T ANEE'US  ESUaTJEAG  usf 
U T>^  SPLvE  e 0R  TrE  CEEEFICIENts  CO  ANT  rs'BN»i///) 

PRINT  u 

Ij  EE-Rma’  (1X..NPT  ALL  IRE  CftEFTLIENTS  E0  genFRATFD  ARE  ylS’E''  :>El£. 

1* I //  ) 

PRINT  yA 

14  EFRmaT  (yX,*AS  A partial  ChECy*  RREN  Nim,  THE  VAyUE  uf  ft  ( E< , N ) IS  G 
1 1 VEN  HY  P. , //// ) 

PRINT  y S 
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U0)».VX,»«(M.S0)«ibX.»D(K,100)*.//l 

t 

t lhl«  P*HT  er  iMfc  PRI^GRAK  Gt^tP*ltS  VF  SrSTFM  3F  SlMUUTASF6Ub 

L k3CAnt)N!>  uStU  IF  SF)LVIFG  ^FH  >Ffc  LetFFUlfcuTS  CO  AND  CSUHF  , 

w 

.rRc»i 

F0IK 

call  HfcS^fcMO  IjERCiNE.x) 
tC  40  »=1,K 
t 1 XX  1 M ) /C 

12»1iO*(lC«Ti) 

T3xPbSj(x(")iO) 

14x(1J..i(l/2.0 

lbxTi»T4 

16XX(M)*A/C 

17xFtSj(T8,C) 

18xBfcSj(T6,l) 

1 5 X (1  ( 1 b«.<  ) 

UCx((L‘-l.C)»X(P)»((A/L)«*4)*7'i.)/(i,0*C> 

F=Tb‘TiC 

IS  » (t<  )x- ( t A-tC  ) • ( A/C  J • HB/X  (P  > ) 

N*1 

?0  IF  (F. !:'■•, •4)  21,52 

21  fc(B,R)xP 

tP  Tt  JO 

22  112XX(M)/X(N) 

113x1. 0/(1, .,-(n2*»?)) 

H4xx  (\).A/C 
T15>xfafciw(ll4,l) 
lltxtfcS-( '1A,0) 

H7x(’ib.Tn-(T12.Tie«^?> 
fc(F,F)i(|.**LC>»IA/(C»./))*ll3*ll7 
OP  30 

30  M\*l 

IF  (F  ,;.c  ,«  ) 20,31 

31  FPIF’  32,-<,  X (“  ) ,F  , Y (►  ) ,F  (►  , 1 1 , t (F  , 10  ) ,p(>' ,b0)  , 3 (“  > loot 

3x  FPffa  (3*,l*,7t3x,Fj>,«)) 

40  Cf■^M^ut■ 

L 

L iPLdlSF  ilk  iHfc  SlSltF  PF  F SIFlL’AFtPlj'  FJI)A^  [FINS 

c (f.L  1-t  C'tFk  iLltMS  CiLB* 

L 

0 iHt  iunkSPlI  .t  FtFL»CFb  TkC  C’FilAFT  SfcCl  PP  "IIm  The  S8LUlIE^ 

C VECie«S,  T«Lb.  The  CFNS1AF3  SfcCl'FS  Y(F)  ofFISFO  IN  STATFNFFT  19 

C APt  FFPLAOt'l  •'T  IhE  S'llIIWF  VtClEHP.  ’hISF  S3L01IPN  vFClbPb 

L APk  lAUb  HELAstLltn  *b  CSlFF. 

C 

CAl'.  •'aTalhIo.t, 100. ICC, O.UFT, loot 


LUi  i 
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c»i-Ci.lat lt'\  a>  uf  r6F»ricifNT  co 

t-ulN’  41 

«i  (////,  1«1  . «C*Lfl.L*’ It  N t»  ’Kt  CutM  irlfcH'  CC*.//) 

F M 1 \ T 4? 

it  ►'‘.‘'4'  (2<1«  . 0M>,  ■ 4.1.T  I fck  J)i  , 12F,  (M'.)**/t),4X|4MY(N),12X,14H>l  JSNI' 

1C.  i."  :►  ) 

F W I M 4 j 

4J  KH»4’  (2''».7HvFCTeos,V>,fH»(v).*/c,b»,*l»,lJX,OM»8FSSELl  = i"*.9FiTr 

1 F “ ^ f F -I  , ) 

F F 1 » T 4 4 

44  F^ff»T  >«>ll>.«T(M4,12x>4XN4C*)ll>i*b:S^CLl*>FX,*T 

lkFF2fcF;..,l()X,.2TtT*t«.//) 

A7£T»Lju,0 
C >t  \,si.i> 

etbStLlsdtbJ(XN4C,l) 

’FFF2kHE4V(M»dtiSfcLl 

4T£T4L:ZTeT4i,4TEK»'2FKF 

FFJM  4bc'F,X(N),Y(N),XX4C.btS«tl.l,1fcR»>ZfchF>,7TflT4L 
4b  f:'ff4’  ( jx  , 1 4 , o<  ix  ,FI?  ,!  ) ) 

■>6  C0NTl\,,fc 

4ls(l‘/C  )••?)•£.  4 
i?»(l.t-((4/(;)»»i) ) ’kC 
43s(4/(C4»2)l*lF|04<L**LC))»27e74L 
F»1M  47,/i.2i’,7i 

47  F0F-4T  (////, b4.jF21»  el3.5,5li3F2k«  E i J , 5 , SX  i 3»<Z  3t  fcl3iS> 

:c*.xi-?2-2< 

FWIM  48, ca 

4t  FPFK4T  (////, b»,3FC0«  tl3,9> 

C C.41.CCL4T  [..ON  l-,t  7FF  T6T41.  4NECII  CL"8fcM 

C 

FtlM  70 

7C  FP,»4T  (////,  1 r.1  , .CAlCM.4T  Hn  tr  Thfc  TbTaI.  4N0D1C  CURHE'x^*  > I 
F&IM  71 

71  Ff'FFA’  ( IX  , • 1 TB7  4L»TeY4t  tixEUIC  CLFKfcM  tIVlOfcO  Hf  T.,;  CeNCl'C-T  I V 1 T 

1 7 ♦ , //  ) 

FOIXT  72 

72  F0FF4T  ( 43X  . oH*  ( M 44/C  1 1 X I 12h.  ( X ( X ) 44/C  ) ,3X  , 14m2  , c«F>I  •4»Ci»iid‘.' ) 
FHJXT  73 

7j  FfH-X’  (6CX,4l«,14X,SF»FfcSbELl») 

FPIXT  74 

74  FPFFxT  (bX,*  »,'iX,4X(K)4illX,4tSLPX«.12X,*XX4C4.1UX.»3E''bFl,14,9X.» 

l7fFFl’.;T«,v>,*lTtT4U4|//) 

MtTXLsO.i 

19  74  .«;,4 

X'.  4C*»(X)*»/. 

Ef  bStL::riFb.J(Xi,4C,l) 

C9l,FX:y(N) 


1,.{  ^ V i 


■ \'-yi 

/“» V i-xJ^aDlL 
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IT6TAu1TH1Al*TF»'>'ITST 

f-RlM  75i‘'.*(M.CSLBK,XKACieES!)bLl.U‘iMlT'“T,P3lAL 
78  fEHMAT  ( J* , I A ,6< JX  , F13  ,8  ) ) 

7?  CPNTl'.yfc 

t 

C CALCELA’ If'N  O-E  LSCA^  Cl/KKEK’  CtNbiTX 

C 

ECIM  t)0 

SO  EPAiyAT  (////,  iHl  , •CAlCCL  A’ JEN  tE  TyJ  L0Cai  Tu’E'^EN'  rENSn»«.//) 
FRJM  al 

«1  FOftyAT  (iX,.iLaCALxL£CAi,  '[vlDEl  t-r  i-'^N '' v C ’ ! V 1 

nx».///> 

810  E^FAC  rtll.wOUT 

ell  EPBE-Al  (H0.0> 

812  REAC 

813  E 0NXAT  (no  . 0 > 

jE  (Rw.Nfc.HCuM  e?0i«c 
8?0  EBJNT  83, «w 

53  EfhyAT  (////, 3x , ORNBi  E?,3) 

EBJM  S4 

54  EORrAT  ( 4Sx  . bhx  ( N ) 40/0  I f X , I2h.  ( X ( a ) • E / C ) , 3 x , 1 i - ( 1 / 0 ) • C A • » ( N > ) 
FBJM  88 

85  EPRmat  (ftOx  . •C*  1 14X  , 9f  *Et  £!>EL?  • ) 

FRJM  efc 

St  EESRAf  (5X,4\»,<;-.  ,«X(A  l^iUX.^tSLRN^.l^Xi^xNWCA.lOX.^^E^bEi.?*,?*.* 
iTERRjLjCtiVX.^lLtCAL*.//) 

UefALtO.n 

EE  «9  N»1 .4 
XNRC*X(N)»kR 
bESSfcL2»BfcbJ(*i'iSC,0) 
tSE=N=x(N) 

lER“lL3CA(1.0/C)»CSusA*X(A)»BEbEEL2 

ltirALilL«tAL*TEh>'lLSC 

FBIM  87,'<,X(M  .CSOPN  ,XKHC,EE5bEL7.TtEMlLPC.lt3CAL 

S7  E ER-AT  (3»,lA,a(3x,ti3,!)) 

S9  CPN’J'.uE 

0''  T£  812 

C 

C CAlCOLATJEN  b'E  Te-E  E£ltKT!AL  r 1 ElE  lei,E  1 I- N 

C Al6A(>  TRfc  "ETAl  bLKEiCt 

V 

9;  EMJM  91 

91  EERmAT  (////,  iHl  , •CA^CI.L  Al  JEN  fcF  Thfc  FElt'TlAL  D I b T R [ 8 y T ; r,  \ . , / / ) 
FFJM  92 

92  EER>-A1  ( IX  , •P0TEMy»t  (R  , 0 ) AElELlRilE  FElt'lTlAL  Ayl■J^r•,  TmE  “Eia,,  S,% 
lE  ACE  •,///) 

920  REAI-  921, Ww 

921  E E»yAi  (H 0 . : ) 

JE  ( R’l  , \fc  . RPyi  ) 930ilC0 


1 
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C 

C 

t 

t 

c 

c 

c 

c 

c 

w 

(; 

t 


933  I M v4 , WH 

94  fPlH^AT  (////,  3*  , 3hRHx  ^5,3) 

F•W^^^  9b 

9i  fPRwAT  ( 4b*  , a«4  (^  ) •R/C  ( « X , I2W,  ( X ( ^ /C  ) , i 9 * , 1 4H«u^^J  I ^r.  Sl“  ?r) 
PRIM  96 

4ft  (43x,.U».13x.l4RCSlti\»  = tbSEi.i:.3*('’H3fcR4R0’L=.3Xi;?R-t3-S^‘' 

1PPTL«  1 
PR ] X T 9b 

47  PPRP4T  (i*,•^•,4*,.X(^)•,:lX,•C,^■^.H^•.l^x,•x^lRC•.lO*.•J£SbP;.?•,8Xl• 
tTf.  hRpf)TL«i9X,«bLPPET|_.,;x,«Pe''tP1|».7/) 
bUPPfc  T|_i(; . 0 

CP  44  , 4 

*NRC»X(6)»HP 

tfciSfcL2*(j6bJ(XiNPC,0) 

CS4R^=V(^) 

TEH-Pf)Tl.iCbUb^»R8PSt-L; 

S4PPtTL’SOMP03L*lERPPpIl. 

PFTPM^:.  (30«SU>'FPTL  ) 

PP14b  98/.  ,X(>.  ),C,SLP\,*XRriFfc5bPu2|3tl'PPiTL.SU4P0TL-P’^TfcN’L 
4e  F PRP4T  (3X,14,7(3X,FJ3,51) 

94  Cl^^U^ufc 
GP  bt  920 


lOC 

SttjRtUTlNb  Btb4FFe(j£RL,\F,7fc':t) 

!3kM  NUMBfcW  - rOOObRCO 

bULfc  - ZF.«f>b  0F  TrF  Ft'SEL  F^NCItN  Ef  7"£  FHbT  Kli<n 
:CtM  N»«lfe  - C3-6Fu-BFb2tFe) 
u4M;i.FUfc  - 36L'0/2?0CI  F(iRTC4\ 

C0PPl.TtR  - CBC-3£00 

CPNTR  ! 09Tl-'N  - jiNFT  F,  pAFisF.  G;LF  7813 

REFEaWCF  CE^P^.T7’Ie^  PtMF-R,  PIS  mviSPN 

8RG4M2*TIaN  - \SL  - \*\,4L  Ffc'fcAFCP  L4HEH4TPHY, 

-4ShUG>EN,  C.r,  iCJPO 
t*bE  • 1 July  1471 

PURPeSfc  - TH  F/L  Th£  first  R fFRES  tF  ,SU9M(X)  FWR  OENEii  WREPE 
P IS  SLPPLIFL,  py  TRfc  lSER.  P ThP  SjHPtiuTIME  CALL 
CIRFNSiCN  X9A0l4),XJAl(3)|X.AJ(2).*9A3(6).XjA4(6)iXJAb(4).ZEP?(l) 
tAlA(XjA0i2.434P<b5b77.',b2C07ell0J.8,6b37279129,li,74ib34439), 

1 l*j4i«3.BJl7Cb9702,7,ci5bP6664e,lJ,17346813b), 

2 I X JA215 . 1356223 .e , 41 7? A41 ) , 

7 (XjA3i6.3HOttl«.S.761C23l,13.rib20c7,16.223464o.l4,4094i4«, 

4 22  ,btJ27?9b)  , 

5 (XjA4*7.b8a3427,il,C64709S,14.372b36/,i7,6lb9660,20.82643XO. 

6 24,3140195), 

7 IXj4s>a.77i4£ 38, i;,33E6C42i15.7001741,i«, 9831834, 22, 2177994, 

“ 25, 4j:3411,?c, 6266143, 31, P117167, 34, 4607813) 

PI»3, 1415926586 
RPLr.«4  , O*J0RU*  JPPC 
GP  7f.  (1,2,3,4,5,6)  jtPL«l 
1 LP  11  1=1. NP 


cooncic" 

0003C101 
0306010? 
00000107 
0006010* 
0006010S 
0006010* 
0 0 0 6 C1C’ 
OOOOOIO" 
00060104 

oooooii'' 
00060111 
0000011? 
00060206 
0006030P 
00000*06 
00060506 
0006C60P 
0 0 0 0 0 7 0 R 
0006060" 
000 6 040" 
OOOOlOflO 
0006110" 
C006l?c - 
0006130" 
00001*0' 
0006150" 
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l>  ( ! , :.T  . « )a'-l  To  iC 

11  CPI^Tl^yfc 

H f ’ I H ^ 

c li 

If  ( I ,0T  . ju,'-'  TO  io 
I ) 

12 

N t M 

J t'  12  Isi.nR 

1 f 1 1 I '•!  1 2 1 ijW  lb  2 0 

2 f » C I ) : X J i 2 ( I ) 

1 i C ' ■^  M f-  ^ c 

X fc  1 L ►>  N 

* C"  1< 

If  ( ! ,til  ,6Ujfl  Tfcl  2 0 
( I ) ax  J*.5  ( I ) 

14  CfN’lMuf 
f-S-  1 l>'4 

1)  LP  15  IaiiV-1 

If  ( I ,‘IT  .6)0')  Tfe)  20 
2f  ><7  I I ) aXJa')  ( 1 ) 

15  CFMINlL 

t Lf  U 1=1.\7) 

l>  I I ,f>T.S))03  ft)  2 0 
I )axoA5(  1 ) 
le  CFMi'.ub 
(-  b 1 1. 0 ^ 

?C  fcfc’»=(f>I/«,J)«(2.0»J(:St..4,0»I-l,0) 
xl=Ht^»»a.O 

x?5xl*.«tl 

4b‘.P(I)!afct4.(HFL[:-l.C)/"l*(i;0*l.C/x?a(«,C*(7,0*M'0Ln-31.'')/2.0 
1 •l,3/-2*<34.0*(83.C*)-l=)LC*x'Lt*?6?,U«iei.C*377V,0)/l5,0 

? ♦l,02x2a(64,o.(e?4i;,C«xSlL«xfLC*x:'Li.-lSJH55,U»H(.-,.D»nPi.n 

a ♦I58574j,(;.x,':;.!).fb7;237.'-l/105,0)))l 

OF  ’£  (11.12,13.14,15.11  )vt)Fi;al 
1C  fci'L 

tU\Cai2)'.  7tSj(X,M 

C4TA(03:,^o27544947be),(ms-.lt5265?t‘.ifc/),(82a.3M83l3l?263f  6), 

1( H 3a-. 9057867427764). (XXs.lGo 3 0196373),  (H5t-, 734853  '5735), 
5(O6a,292l?672487fc-2).()'7a.,65x5ci7ta,7i^.8),,KH,,4453H0i«051F»»), 
'(SQa, ?d 2784494766), (Sl«. 21 69 554774 376), (5?a. 700468 2514763), 

4(40X4. 532  ' 42093762).  (41a4. 22177:411661), (42i5,2443314672E-1), 

'7  (50=, 4488459489663).  (Fla, 75725C4657962), 

‘(FCa- 1,23  3 944555161), (l.la-2, 7767921059), (F2a-4, 95173991266-2), 
7(rii,4i00854523'-2),(f:64,).  (014,722  3 6648621). 

» ( :0= . 17496x/o23963 ) . 

4(4xr!i. 9838/27495977). ("61  = -. 1142532572177), (H82=, 4094621362565). 
-*(=5  3 i-,r6!S6Jil985663).(784s, 5757541403561),  (885  = -, 279  ■J4475519F-1), 
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riS!j7s,iVt,l/4!>4«9  1to),t!>'l!,16711t7Jlb4fc6),(5!>i?«. 60777  ?6«i!4  7FJ), 
F(»-b0*,62«J6HS66  71F3),(oFl:,V730C0946£et2)> 
f(ri,0».2U«b47«i31F3t.(i;Cl=,46':i7l2?6i<:t2). 

n(  440  = 3,5>4blbVVsl7iF2),(»Ai:S,b':44fc4J02ltl).t*A2i6,b2P3Cl8«4«t)f-l). 
►•(CC.C44,4B?2J4H24t£l),aCltV.774ec6fc76  4),(rc?*l,  7 728b  7 7 1 4 jF  . ) ) 
t :X»* 

liti  Ti^  6 J 07  ^6  7 » Qti  73  e 

6 17 (C->  )1,1,2 

1 F j I ( ( ! nt)» 7 > ‘riRB  1 »r ♦Ki> ) »r ♦9  4 ) »r  t p« ( i (p^r3 )*o*h? ) ) •l*rc 

faEij  =P/ ( ( ( ^*02  ) *r  ♦EC  ) 4 FfOFN 

2 If  IC.r.T  ,(,)  G0  7 7 9 

A = APi(x  ) »■  Os'  /C 

f = i(*2«c^4i).o*A0)/((r^Fi>*i:^90) 

•(C0i>(A)^(P.C)^bIMA)*(P-C))/bCR7(4) 

F b T L b f. 

7 If (C-f  >11.11,21 

11  F«(  ( ( ( ( (RR7»U^Rbt  )«n»Ck5).iJ*RC4).C«Bh3)«4,.BH2)*r^‘'Rl)»0^f'PC 
fctSijAX.b/(((D^bS2)#n*fbl)*i,»S'C)  ? kfcTu«N 
21  ifir.GT.Gi  G0  rr  9 
AsAHSj(X)  * Dxf/r 

f:((4A2.3»a41>»C*aaO)/I(I;»o‘-H*''»PcC) 

C=l(CC2’J*kCl>«r»CCO)/(A»((C^''il)*fLro>) 

Ar(  C7b( A ) ♦ ( J-P  ) ‘SIS ( A ) ♦ 1 C ♦f  ) )/ECf  ’ ( 4 ) 3 Irtx,.'  ,U)4:-4 
b b b M * 4 

(.  b 7 0 b ‘•I 

c bS!\’  <jl,'« 

8l  f 7b“AT ( //ijx.t B9tc  IN  Obb,.  N =*15) 

G'-'  7t  ;J0 
9 FKl*  ■’9;.x 

91  f fbPA7  ( //l,x.tBO(:K  IK  obS.i  A9GLbEN7  X T^P  L4»3b.  X : ♦ b 1 7 . I C ) 

ir:  bbb^=i.b3:c 
tNC 

bUbReU7lNF  bA7Al.C(A.X,N9,NV,irt,LE7,N4C7) 

LlbbNblC-.N  4C.AC7.NACf).X(NAC7,NAC’> 

I f ( I L -I ) 1,2.1 
1 Lb  .5  Ul,-.b 
LP  4 

4 X ( 1 , , ) jc , a 

3 X ( 1 , j ) 5i , a ■ 

NVsNb 

A Lf ' = 1 . 0 
NWliNb.l 
i"  5 A5l,  -.Rl 
I '•  1 = b • 1 

b : =c.  a 

LP  6 U«,NR 

X : ACbb ( 4 ( I ,<  1 ) 

Il(7-P;vii7)  s.6.7 


BbSj 

17 

3EbJ 

1< 

abSj 

I"* 

absj 

16 

or  S J 

17 

6FSJ 

1» 

of-Sj 

19 

eFSJ 

2"' 

Bf  sj 

21 

BbSj 

22 

ofbJ 

27 

bF  S J 

24 

fif  Sj 

2b 

8PSJ 

2* 

tJfSv 

27 

BFS  J 

2R 

bb  b^ 

29 

8 = Sj 

30 

BbSj 

31 

BbSj 

32 

BFSJ 

37 

bb  Sw 

34 

Bf  So 

36 

Bf  Sj 

3* 

OFSJ 

37 

BFSJ 

3* 

SESj 

39 

8ESJ 

40 

BFSj 

41 

oES  J 

4? 

B f S J 

4 X 

BFSJ 

4 4 

BESJ 

40 

BESJ 

46 
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OOfc 
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MVRT*^ 

IPmi 
CR^T INot 
IF  (PIVflT  ) 

cei>o,o 

H t T L F«  N 

irut-s.K)  xJ»ll*10 


tf.Vto 


i;  tf>  i<  j»K.NM 
i«*( IPH. J) 

«(1PR>j)iA(K,J) 

U A(N,.)l/ 

tF  13  J»l.NV 
i*X( IPH. J) 

X< |Ph.j)«X(K,J) 

13  XCHi.Fti! 

LtTfUtT 

11  tFl«LbT»A('F.'0 
p 1 V P 1 « 1 , C / * I X • > 

10  14  J»1W1.NH 
4(K,,)sA(K,J)«HIVf.T 
DP  14  i«:mi.NW 

14  A ( I )>A  ( I , J )>A(  1 .X  )«A  (K  ,.  ) 

DP  5 

IMX(K.J))  is.i.l* 
lb  X(x,,)»X(x,jj»hIvpt 
DR  It  l»I«l.N« 

It  X(  1 ,»)!«(  I . J)-A(  1 ,X)*X IX ,,  ) 
s CRNTINjt 

If  < t<1m,  ) 17,i,l7 

17  LtTtDtT*A(NHiMHt 
FIVPTxi,0/A<NW|XS> 

LP  le  J»llNV 

X ( NR  . J ) »X(  N‘<»  J1  *P  1 V0I 

DP  It  K«l.NRl 

I«NR>X 

lUPxO . 0 

D6  IV  l«I.>H1 

IV  SI'X«hli“*A(l.l*l>»X(L*li.) 
le  X(1 t.)sX< I , J)>bix 

tNC 


016 

017 

016 

019 

020 
0?1 
022 

023 

024 
026 
026 

027 

028 

029 

030 

031 

032 

033 

034 
03b 

036 

037 
036 

039 

040 

041 

042 
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044 
C4b 

046 

047 
046 

049 
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051 

052 

053 

054 
05b 
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Appendix  C 


COMPUTER  PROGRAM  FOR  COPLANAR  CONCENTRIC  CIRCULAR 
ELECTRODES  WITH  UNEQUAL  POLARIZATION  PARAMETERS 
UNDER  THIN-LAYER  ELECTROLYTE 

l-f.tGt'AM  yNtOTMJK 

THIS  ewBiiWAM  Cb^-fl,TES 

(1)  CLhRfcNI  LISTHIBLlUK 

12)  P?1EMl»u  ClblRItintN 
(3)  AMiCIC  CuRt-EKT 

EFR  CRNCfcNTRlC  riPCDL»R  fcLtCTRErfcS  C6»/b''t!:  RY  TmIN.l*Y£;<  FlECTRFLYTF 
E0R  IHfc  C*bE  where  ANFtIC  ARC  tATt-tEIC  waFsER  PfeLARIZATItF 
rahareterS  are  net  fclal, 

AsRALlub  0)  AN0CE 
CsRAtlyb  «E  CATweCfc 
RjtIbYANtt  AlWNF.  RADlLb 

ELEMiTRiTMlCuNEbb  ©E  ElECYRFlyIE  LAYER 
RR*R/C 

EA1AN0U1C  *AGNEW  P8LARJZATIFN  EARAREIEr 
tC*CATHfiDlC  wAUKER  PeLARIZATlEN  RAWARETEr 
XIRIsMte  ZE«0  0F  HESbEL  FLNCMfcN  1*  ERufc"  J0PD 
P(R,0>»lMTtR)  ACIAL  PEYEMIAL  A1.FNG  YRE  REYal  SjREAcE 
E(R,O)tP0YENTL«ELECYErt  PFTESTiAL  ALERQ  The  PETAL  SURFACE 
E <R|  0 >iC0MbTANT.E  (R,  0 1 t pITR  YRE  CENEYaNY  caNCEllInE.  3uY  IN  The 
FINAL  EXPRbSbl0N,  SP  CEASYaM  IN  EFFECT  CAM  BE  SET  E3JAU  T0  AERP, 

THE  CBEFFlUfcNTS  CO  ANL  CEUPN  ARE  CEFInEU  IN  TRE  F0LlPwINF,  ESUATITn 

K 

F(H,CI  ! C ♦ SLR  (C  . (X(N)*F/ri) 

U Nsl  N 0 

lTtYAL»TaTAL  anFLIC  Cl-RRtN'  DIVILEL  BY  7h(.  r0MDUCTlVITy 
lLtCALsL0CAL  CuRRFnT  CENSITY  ElVlTtC  EY  T H£  CBYPUCTIVITY 

real  LA.lCi ITBTAl , ILECAl 

tlRENblEN  »(1P0)|Y(10C, IOC), 0(100,100) 

FI»3,l41b92T>!>36 

I "*100 

1 REaC  10  , A,L,o0EN0Yp,l a ,lC  , t a ,F L 
I 1 0 * PRRA  T ( ;F  , C ) 

I PRINT  i1,A,C,hlENGTh,LA,lC,EA,EC 

j H FPRhaT  (ihi,sx,2ra,  F 1 C .b  ,b»  ,5RC*  F 1 0 , t> , 9 Y , Rwd.ENGTH!  E 1 i, . 9 , 9 » , 'H^ 

I 1»*  Fia,E,bx, jmlC*  F10.9,9Y .3NRA*  E 1 0 . 9 , 9 * , 3 wEC « Flo. 9,////) 

1 PRINT  ii 

‘ li  FFWRAT  (i*,«gF_nFRATI6N  EF  TPE  bYRYER  FE  blMULYANEHUb  FOUAYIFNR  USF 

U YF  SblvE  FaH  YpE  rEEFFlClEN’b  CO  ANC  Cbl'BN*,//’/) 

I PRINT 

I )J  FTPRaT  (j«,.\0T  all  TPt  CFtFFICIEN'S  90  OfNFRATED  A»E  LlbY£r  "EuE- 

I !•,//» 

1 PR INT  jA 

lA  FFEPA’  (i«,»Ab  A Partial  ChEC«.  wren  n*R,  the  YAl^E  B''  HIH.M  in  ’■ 

I llVEN  HY  P«,////) 


I 

I 
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t-SIM  lb 

no)*»vx,«f’(M,boi*.ax,»fc(X',ico)».-'/) 

c 

c IHIS  PA*<r  af  ImE  PHOGRAX  CtKfcRATtS  ThP  SYSTEM  3F  SIMULTAStEub 

c tEUAiioNS  usfcu  lA  sn^viAa  ffr  't-fc  tEtPr  icients  cq  ano 
c 

.ppr=i 

\P  IK 

C‘LL  HfcS^twR  (jFRr.NE.X) 

LP  AL  >A  = 1.X 

XRbC«X(M).HLEMaTt-/C 

ilAv-xRac«(txp(xYfcC)-txP(*xMtcn/i.o 

CEiF'XHHCi(tXP(xY6C)«EXP(-XKPCn/x.O 

lAAKXMaC»SlNHXMPC/CPSE«YriC 

n«xiM)/c 

72  = l,0»(l.C«Tl.TAAF,xHbC) 

73«EtSj(x(M)  ,(') 

14*(13**21/2.0 

Tb=Ti.T« 

1 6tX (R ) • A/C 
17iRESj(Tb, J) 

1B»FESj(T6,1) 

1<>s(77«  «2  )♦(!«•  «i) 

nC=<(cA-CC)*X(Y).((A/C)»»ii)»Tt.TAA>XYBC)/(?,0»C) 

P=Tb*T lO 

19  Y(Y)  = -(l,0/CebHXMC)»(EA-FC)«(A/C)*(lF/xc  )) 

N = 1 

?c  IF  (A.b:;,R) 

?1  b<B,MsP 
()0  Tt  JU 

22  11bxX(M)/X(N) 

713  = 1.0/(1, j-(71i..2)  ) 

7 1A  = X ( N ) .A/C 
7 lb=BtS,<(  7lA,  1 ) 

716=fa6Sj(7lA,0) 

7 l/= ( 7ib. 7 7 )- < 7 1 2.T16* 1 e ) 

XAbC*X(A).tjLbNu7b/C 

bIAbXNdCi(bXP(xAbC)-bXP(-XNbC))/x.C 

b(Y,A)s(LA*l-C)»(X/<C».<))»(l.C/CeSt-xbPc)»7l3.7l7.!>lMHXN8C 

OH  7£  jc 

30  A=N»1 

IF  (A.lE.X)  iU.3l 

31  PPIA7  3b , 1, X (P ) ,F  , Y ( Y ) , E (Y  , 1 ) ,b ( Y . 10  ) , a ( Y, 50 ) > 3( M*  100  1 

3b  FPhYA7  ( 3X , 1 4 , / ( 3X ,bl3 .«  ) ) 

40  C0N71\ub 

c 

C iFLUIo;.  •■)►  I«b  SYb’tY  EF  4 b I Y C C 7 A A b F Ub  FOUATIRNb 

c FFb  7Ht  C'-bFFICIcFTC  CbLbA 

c 


i 
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L Tht  bU«t<aunNfc  RfcPLACFS  Tt-t  C*'^ET^^T  VtCTfRS  wPH  THE  SHLUTie^ 

C VECTtHs,  TMoS,  THE  CEN'TaM  vECTF*<S  » ( I DEMNED  I ST*TEkem  19 

C AWE  HEplaCEO  HY  the  SELLMW^  VECTEHb.  ThFSE  S3LUTIBIY  vECTSRb 

C AWt  L^TbH  HElABELLtn  Ai  CSUFN. 

C 

call  HATAUi<h,y.lOO.lCC.O,UET.lCC) 

C 

c calcllatjun  0^  The  cc 

c 

F « 1 N T A 1 

41  FPHhat  (////, iHl  ..CAlClLATIEN  £F  the  CfiFFFICIEyi  Co».//) 

FF!\T  4b 

4ii  FPHhaT  ( iJ8X  . bHiELUT  1 EN  ( 4 3Y  I 12H«  t » ( M • A / C ) , ?X  i 1 0 Mb  [ HH  ( X < M • . 7 K 1 1 3h 
1 V ( \ ) »S 1 AwX\HC I J» , IAHKL NA ] A u SLF  EF) 

FR1\T  4l 

41  FPHHAT  (4[hx.7wvFCTERS.«Xif"»(')**''Cie»,»l*,l7X»10H8LE9GFH/C)»?X,9H 

14HESbFLl»'/»t’HTERHZEFt*J 

F R I \ Y 4 A 

44  Ff-hMjY  (SX,«\»,9*,.*(A  )«,ll»,»»(A)*il?Xi4>N4C*.ll*t»b?SSEl  l»i9Xi*S 
llNHXAnC’>7Xi»TtHF7EHE*<lO)'<*/’CT4L’>//l 

ZTEtal  =3,0 

CP  4t  \=1,H 

X'.AC  = *(')*A/C 

eF  bREul  = c(Fi>J(  XNAC  , 1 ) 

X bcO  = T ( '. ) •3le'»u  T f / C 

bINHX-.„i«(c*HU'cri-t»'-(.Y\FC.))/i.C 

1EHFZFHl-iY('j)»oFiRELl*bIbF*ABT 
Z T E T A L ! Z T p T A L ♦ T E R A' Z E K F 

FR[AT  4b,f'i|X<N).T(M,XNAC.obb':tLl>SlNHxNor,TEHTZERB.ZT0TAL 

45  f PH  hat  (3X,U,7(Jx,fiT,5)) 

46  CP^YI^L^: 

Zl= ( ( A/0  ) ..2  ) .EA 

z2s(1.0*(<*ZC>»»il)*EC 
Z3i(A/(C»»2))M2,0»(LA'LC))*ZYtTAL 
FRlbT  il,nu2,2i 

47  (FHh*t  (////,  5x  , ZHZ  1 : 1 1 3 . 5 I b>  1 3F  zit  = E l3  , b , 5X  i 3hZ3*  E13ib) 
C0»-Zl-/2-z3 

FH IF  Y 4r,C0 

4c  FPHH4T  (////, bx,ZHCO=  El3,51 

c 

c calcclayio\  t'F  tee  Total  4^ECIc  cirhekt 
c 

FW1\Y  70 

7;  fthhai  (////,  1m]  , =0X00  cl  a Y 1 EN  EF  TFF  TpTAL  ASaDlC  CURREbiT  • , // ) 

FRIA I 71 

’I  Ff'HHAT  ( IX  , • I TuTALsTfTAL  ANECIL  CLHREKT  LlVIUfcO  THE  CEKCUtTIVlT 
I Y • 1 //  ) 

F R 1 A ' 7 ? 

fplhat  ( 45X  . B«*  ( a )•  A/C  .6X  , 1cm.  ( X ( a ) • A/r  ) ,2X  , 13ms  1 SH(  X (N  ) • . 7X  , llM? 


I 
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f-RlM  7i 

1i  ^R^^'AT  (ftO*  . . 17*  , iCHELFNGTl-/C  ) . 17K*SISMX^RC*8(:SSE11  ) 

M<  JM  7< 

7 4 fe^.^AT  (;)*,•^•,7»,«)^(^)•,ll*,•CtLF^•.ll*,•*^;AC•.llX.••lt'^b^Ll•.9».• 

lTeTAl.»0,3 
to  79  N=l,n 
*^AC»  * ( ^ ) • A/U 
fatijSfcLIxbFbJfXNiL.l) 

CSLHN*v(M 

*^bC«*(^)•BLfc^bTI■/C 

S INRXNbCi (fc*F( XXfcC ) -tXP(-XNEC))/i.t 
^tHFlT0Tl^,J.Pl.A•CSLF^•bEbSfci.l•iI^R*^Ht 
1 ItTAUi  I TP)TAl»TFFF  I Tt  T 

FRJM  7b,-'i,XtM,tStH*,XXAC.FE:FbFHiSlKMXN-C,lE8“IlBl.ITPlAL 
7b  E f RPAT  (JX,M,7tiX,El7,*)) 

79  CONTINUE 
C 

C CAlCLLaTIKN  0F  T|-f  LCCAL  CURRFM  rfcXblTV 

c 

FRlM  dC 

PC  PFRFAT  (////, 1h1  ..CAlCUAT  ion  tF  It-F  tec»l,  CuH'lfcNl  Dt  JSllT.,//) 

F F 1 N 1 b 1 

PI  FOrfAT  (i*,«lLaCAL*CCr*L  CuRNFM  CtNiliV  •divided  «v  T^ae  CONPUCTIVI 

ITT.,///) 

bic  feat  Hii,«(;uT 

8U  E PRfaT  (HO  . 0 ) 

612  R6AC  fiU,RK 
ttlj  EGRFAT  (flC-OI 

IF  ( RR . NE . RCuT ) t ?0 I 9C 
6/0  friM  eb,FR 

pj  Ff^RFAT  (////,  3X  . JfRRa  E;,7) 

print  6A 

64  FFRfAT  ( 4b* . 8HX (N ) «R/C »<X . 1/R.  ( * ( M • R /C ) . 2 X , 1 0”$ 1 Nh ( X ( N ) * , 7 * , 1 6h ( 

ll/O^CbLdN**!**)  ) 

RR  U 1 ob 

Pb  fORFAT  ( 60*  ■ .0.  . 17X  , ICfElFNGTr/C  ) . iX  , 17R«RInhX'AHC*3FSSEL2  ) 

FRJNT  Sfi 

66  F(^RFAT  (b*.*N».9*,.X(N)».llX.«CSLFN».12*i.XNRC:*.lUX.*lESbEu7*i9*i* 
lblNF>NdC»i7X,*iTtRFlLEt..9*,.tLfC*L..//) 
iLtCACsO.O 
C^  P9  >=1,R 
»>  Rr«X  (\).RW 
rE  bSt ^2  = d6b J ( xnRC  . 0 > 

CStfcN=v(N) 

XNbC'X  (N  ) *BLENIjTF/C 
bINR»NHCl(6*F(*Ncr)-tXF(iXNFC))//.C 
TFR";i-c)(;«(l.C7C)*CbUbT**(N).S:NF*NeC.FtbbKL3 
IltOAus!Ll-CAi..TtHFILcC 

RRJ\T  e7,N,*(N).C'CPN.«NRCitt'b>-L?>SlNR*NRC.TfcR“IU^C,It'’CA.. 
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C 

L 

t 

C 


07  F0M'*T  ( j>,  14, /(JX.EU  .5  ) ) 
&P  Tt  al2 


CALCULaI  I'-)N  7>.fc  PtTfcMlAL  riSlKlfaLTlEN 

ALtKC  Tl-fc  kETAl  SLHEaCE 


9C  Pf>l^T  91 

91  fPhM*T  (/■///, Irtl  .•CAlCELATJEN  tf  TOE  E0lE''TlAL  0 I ST  R J 9|jT  [ BN  • , / / ) 
PRINT  92 

92  fPNH*T  ( IX  , •PBTENTi.st  ( R , 0 ) sturtTRf  CE  P0TtNTI4L  AGEING  ThE  *'E’Ai.  SL° 
IE  ACF • , ///  1 

92C  REAH  «21,HH 
921  f BRPAT  (HO.  3) 

IF  (RR.te.hCuT)  930. ICC 
930  PRINT  94, RK 

94  PPRHAT  (////, 3X,3PRRx  P«,3) 

PR'iNT  95 

?b  PPRHAT  (43X.12M.  ( X(N  )»R/C  ) .2X  , UPCEEM  X (►■).,  7X.  i3hCS  JBN»l'ESiPL2 . 3 
1X,14PRuNMnG  Su“  PP  1 
PRINT  96 

96  PPRmaT  (44x,«O».17x,lCPELENCTP/C),3X,loP»C0SMX9yC*.6X,9MTERMP»TL«, 
16X,liH-Co-bUMPaTL»J 

PRINT  97 

97  P0RMAT  (SX,»N«,9X,«X(N).,llX,*tSLBN«.lOX»»BESE£L2*,9x,«C0SHXNnC».A 
IX  , •TtRMPaTi.*  , NX  .•SO*'PETL»i<TAi»PPTENTi»,//) 

SEPPeT^.so.C 
L0  9«  N»1.N 
XNRPsX ( N ) (RH 

bESSEL2oBEi>J(XNRC,0) 

CSEPN»Y  < N) 

XNbC»*(N)»BLENUTP/C 
C0bPXNHC»(tXP(XNEr)‘EXP(-XNRC) )/< .C 
TERPFWTEiCSUBN»RESSEL2»C6SPXNP(, 

bUPPfcTi_sauMP0TL»T£HMRE  IL 
PPTENT|_s.(tJ*Su''PPTL) 

PRINT  9e,N»X(N)  .CSUBN,BESStU2.tEEPXNbC,  TfcRHPST,..Sui4P0TL,RPTPNTl. 

9e  PPPPAT  (3X , 14, 7( 3X ,E13  ,«  ) ) 

99  CPNTlNuE 

OP  re  920 
lOC  ENU 


SUBROUTINES  BESZERO  (JORO,  NO,  ZERO),  FUNCTION  BESJ  IX,  N),  AND  MATALG  (A,  X,  NR,  NV, 
IDO,  DET,  NACT)  ARE  GIVEN  IN  APPENDIX  B. 
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